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Abstract. Current status of (the lack of) understanding Blazhko effect is reviewed. We focus mostly on the various 
components of the failure of the models and touch upon the observational issues only at a degree needed for the theoretical 
background. Attention is to be paid to models based on radial mode resonances, since they seem to be not fully explored yet, 
especially if we consider possible non-standard effects (e.g., heavy element enhancement). To aid further modeling efforts, 
we stress the need for accurate time-series spectral line analysis to reveal any possible non-radial component(s) and thereby 
let to include (or exclude) non-radial modes in explaining the Blazhko phenomenon. 
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INTRODUCTION 

Efforts in understanding the Blazhko effect (periodic 
amplitude- and phase-modulation of the pulsation of 
RR Lyrae stars - iflrjll . (Hst) date back nearly to the 
discovery of their variability itself 11511 . Unfortunately, 
the situation has not changed much since then. We are 
stunned by the overall simplicity of the phenomenon and 
the lack of theory that is able to account for the basic ob- 
served properties. Although several observational facts 
indicate that the theory of horizontal branch evolution 
and RR Lyrae pulsation is basically correct, and it is un- 
likely that the future understanding of the Blazhko phe- 
nomenon will involve dramatic changes in these models, 
it is still highly unacceptable that we do not fully under- 
stand this basic class of variables. We agree with 112711 : 

"Explaining periodic, or nearly periodic, long term 
modulation of the light curve seen in many RR Lyraes 
- known as Blazhko effect - constitutes perhaps the 
greatest challenge to stellar pulsation theory." 

In this review we briefly summarize the various efforts 
that have taken place since similar reviews of the topic 
(H, see also S and (H]). We attempt to give a 
simple discussion of the reasons leading to the failure 
of each of the current models. As a matter of curiosity 
we also present a toy model that, for the first time, shows 
a behavior similar to what we observe in Blazhko stars. 
We conclude with suggested directions of research and 
expectations from satellite data. 



OBSERVATIONAL STATUS - A QUICK 
GLANCE 

RR Lyrae stars 

There have been several thorough works done dur- 
ing the past seven years that have substantially con- 
tributed to establishing a more complete picture on the 
observed properties of Blazhko-type pulsation. In addi- 
tion to the analyses of the large pho tometric surveys of 
MACHO, OGLE and AS AS II, H, Hi, there are ob- 
servational projects focusing on the particular problem 
of Blazhko modulation. Johanna Jurcsik and coworkers 
at the Konkoly Observatory have been monitoring these 
objects photometrically by a dedicated 60-cm telescope 
since 2004. "The Blazhko Project" led by Katrien Kolen- 
berg at the University of Vienna studies specific targets 
both by photometric and by spectroscopic means. Fur- 
ther photometric studies are made at the Michigan State 
University by Horace Smith and his students. Very im- 
portantly, there are spectroscopic investigations carried 
out by Merieme Chadid, Katrien Kolenberg and cowork- 
ers aiming at to measure magnetic fields in RR Lyrae 
stars. We refer to the following home pages on the works 
done by these groups Q 

The main photometric properties of the Blazhko stars 
are very simple (as compared to some familiar classes of 
multimode variables, e.g., to those of the 8 Scuti stars). 
The frequency spectra of the light curves are (usually 
strongly) dominated by the pattern of 



Vbl, {kv Q }, {kVo±v BL } 



1 http://www.konkoly.hu/24/ 
http://www.univie.ac.at/tops/blazhko/ 
http://www.pa.msu.edu/astro/smith/ 



where Vo is the main pulsational frequency (that of 
the fundamental [FU] or the first overtone [FO] mode), 
Vbl is the modulation (the Blazhko) frequency and k — 
1,2,..., usually goes to fairly high harmonics in the case 
of the FU (RRab) variables. The amplitudes of the modu- 
lation components corresponding to {±Vbl} are usually 
asymmetric: 

A(kv - v BL ) ^A(kv + v BL ) 

The highly asymmetric cases often go to the extreme, 
with one side under the detection limit - e.g., [57]. These 
cases lended further support to the non-radial mode exci- 
tation concept. At the same time, the explanation of this 
asymmetric amplitude pattern (even if it is only a mild 
one) is the main stumbling block of the current models 
invoking non-radial pulsations. 

In a non-negligible number of cases the frequency pat- 
tern is more complicated, starting from the appearance of 
higher-order modulation side lobes (see [35], [40], [29]) 
and ending with more complex structures around the pul- 
sation components (e.g., [87], [53], (HH) or even more 
intriguingly, with the appearance of very low-amplitude 
structures with the same modulation frequency but cen- 
tered at 'no-clue' frequencies (see Chadid et al., these 
Proceedings). All these add to the flavor to the Blazhko 
phenomenon and could be instrumental at some stage 
in searching for the physical model, but the more gen- 
eral pattern described above constitutes the observational 
framework the models should deal with. 



TABLE 1. Statistics of Blazhko stars. 



Type/FielcT. 


LMC 


SMC 


Gal. Bulge 


Gal. Field]. 


RRab 


12% 


10% 


25% 


>5% 


RRc 


10% 


10% 


12% 


>7% 



* S ources: |701(SMC), |1](LMC), |46](SMC), |53](LMC), 

1 19](Bulge), |77](Gal. Field), (5j (compilation) 

' Note: See comments in text for the expected increase of the 

incidence rates for Galactic RR Lyrae stars and for the estimation of 

the rates in the case of multiple sources for the statistics. 



The relative number of Blazhko stars varies in differ- 
ent stellar populations. Table 1 attempts to give a current 
census. When multiple sources are available, we present 
the values considered to be the 'most reliable'. We cau- 
tion that considerable changes can be expected in cer- 
tain populations due to the advance of measurement and 
data analysis techniques. This is especially true for the 
Galactic field population where both the ground-based 
survey of Jurcsik et al. (these Proceedings) and analysis 
of a handful of CoRoT RR Lyrae stars by Szabo et al. 
(these Proceedings) suggest that the incidence rate could 
be (even) higher than 50%. Obviously this is an impor- 
tant issue from the point of view of model builders and it 
will be interesting to see the distribution of the modula- 
tion amplitudes; whether they follow a continuous long 



tail to the micro-modulation regime or they show a bi- 
modality with peaks at the 'classical', high-amplitude 
Blazhko modulation and in the micro-modulation range. 



Cepheids 

Before 2004, except for V473 Lyrae (HR 7308, see 
Hill for the analysis of the radial velocity and If37ll of 
the Hipparcos data) no Cepheid variable was known with 
modulated light variation, similar to those of the Blazhko 
stars. The analysis of the OGLE database of the LMC by 
Moskalik et al. f50Tl has led to the discovery of a substan- 
tial number of Blazhko-type Cepheids. This is a very sig- 
nificant result, since then we need to seek for models of 
similarly modulated pulsations in two classes of stars that 
are dissimilar to a great extent. In Table 2 we summarize 
the current statistics on these stars, based on two avail- 
able sources 17 111 and 115 211 . The two studies yield basi- 
cally consonant statistics, in spite of the different datasets 
the works are based on ( 17111 use mostly OGLE III, that 
is more extended and more accurate than those of OGLE 
II of 115 211 - however, 17111 do not give any details of the 
analysis but only the final statistics). The systematically 
lower rates in f52ll can be attributed to the less extensive 
data used in their analysis. 



TABLE 2. Statistics of modulated Cepheid variables in 
the LMC. 



Type/Ref 


FU 


FO 


SO 


FU/FO 


FO/SO 


EU 


4% 


28% 


28% 


18% 


36% 


[52] 


0% 


8% 


12% 


13% 


> 19% 



We see that modulations are observable in all types 
of Cepheids in quite high rates. RR Lyrae stars are ap- 
parently different in this respect. In the case of double- 
mode (FU/FO) RR Lyrae stars the Blazhko behavior has 
been reported only in one object 115 311 . Interestingly, the 
phenomenon seems to be quite common among double- 
mode Cepheids, especially among the FO/SO ones. It is 
important to note that the average of the modulation level 
A mo( }/Ap U i s of Cepheids is only about 5%, i.e., substan- 
tially lower than those of the RR Lyrae stars - e.g., ll77ll . 
Except for the FO/SO stars, which have modulation peri- 
ods greater than > 700 days, all other types of Cepheids 
span a similar range of modulation periods as those of the 
RR Lyrae stars (compare HI and l52ll ). A very interesting 
feature of the modulated FO/SO variables is the oppo- 
site phase by which the amplitudes of the FO, SO modes 
change (in the modulation interpretation of the observed 
frequency pattern). This is certainly at odds with models 
based on geometric effects in explaining amplitude mod- 
ulation (i.e., the oblique rotator and the 1:1 resonance 
models). Some dynamical mode energy exchange model 



might be more appropriate (as mentioned in lHH), but 
there is no idea on the specific physical and mathemati- 
cal details of this interaction. In any case, these new dis- 
coveries may perhaps open new ways in searching for 
additional ideas/models of modulated stellar pulsations. 

CURRENT MODELS - WHY DO THEY 
FAIL? 

In this section we briefly describe various ideas to ex- 
plain modulated stellar pulsations and then discuss in 
some detail the three, currently most often cited mod- 
els/ideas. We stress that there are very important differ- 
ences in the level of elaboration among these models. 

In terms of elaboration the non-radial resonant rota- 
tor/pulsator (NRRP) model stands out, since it is based 
on the detailed analysis of the nonlinear interaction 
among radial and non-radial modes (see ll83ll and JUl]). 
We recall that nonlinearity is a basic attribute of the 
Blazhko phenomenon, since it is observable throughout 
many decades, i.e., well above the dynamical (or even 
thermal) time scales. As far as the methodology is con- 
cerned, the drawback (and, because of its generality, as 
well as the power) of the NRRP model is that it is based 
on the amplitude equation (AE) formalism (see l23l) . i24ll 
and |01 JE and, as a result, no strong specific predictions 
can be made for the observed behavior of particular stars. 
Nevertheless, the predictive power of the NRRP model is 
not worse than any other available models/ideas. 

The magnetic oblique rotator/pulsator (MORP) model 
lacks the complexity due to nonlinearity and does not 
examine the stability of the presumed non-aligned mag- 
netic, rotating and pulsating configuration. Instead, it 
postulates the stability (perhaps based on the observa- 
tions of spotted stars and rapidly oscillating Ap stars) and 
gives observed properties in a pure kinematic/linear pul- 
sation framework. 

The 2:1 resonance model of [48] is also based on the 
AEs and assumes radial mode interaction only. This is 
the first (and so far the only) model that represents the 
modulation by dynamical mode interaction. It is reason- 
ably elaborated (further works were probably discour- 
aged by the lack of hydrodynamically confirmed numer- 
ical computations). In this review we construct hydrody- 
namical models that may suggest some level of viability 
of this idea. 

There are also ideas that consider some sort of en- 
ergy transfer hypothesis among some pulsation modes 
to aid amplitude modulation. The hypothesis of 12711 as- 



2 Due to the extreme level of complexity of 3D time-dependent long- 
term hydrodynamical simulations, the only method currently available 
to us to study nonradial nonlinear pulsations is to use AEs. 



sumes that there is an energy sink originating from many 
low-amplitude nonradial modes. The interaction leads to 
lowering the amplitude of the radial modes, that even- 
tually results in a modulation of the pulsation (through 
some, so far unspecified interaction between these two 
classes of modes). The authors put forward the slight ef- 
fect of the overall lower amplitude of the Blazhko stars 
in the Galactic Bulge as compared to that of the single- 
period stars (see also the exhibition of the same effect 
for the Galactic field and M3 variables in ll79ll and for 
the LMC in (!]). However, a simple resonant mode in- 
teraction can also lead to amplitude decrease (see [8J]) 
and there is no guarantee (or any obvious physical rea- 
son) that such an interaction will lead to modulated pul- 
sation. The argument by [52] for a possible dynamical 
interaction in explaining the amplitude modulations of 
the FO/SO double-mode Cepheids is more convincing, 
since the data can indeed be interpreted as the result of 
such an interaction. However, the authors stay vague in 
specifying the type of interaction that may take place in 
these objects. 

The idea of multimode non-resonant interaction also 
came up in [42] as a possible source of amplitude modu- 
lation (for the analysis of the three-mode non-resonant 
interaction see also [31] and [86]). On general physi- 
cal ground, it is easy to show that more than two modes 
are necessary to expect non-stationary (i.e., non-constant 
amplitude) asymptotic states Although an extensive pa- 
rameter search may lead to finding modulated solutions, 
it is nearly impossible to construct one in which one 
mode strongly dominates the others (as observed in real 
stars). 

In a current idea presented in lf75ll we find a loose 
description of the highly complex scenario of mag- 
netic/convective interaction with radial pulsation as a 
possible cause of modulation. Although we cannot ex- 
clude that such a complicated system may also lead to 
the modulation of radial pulsation, there are some mis- 
interpretations present in that paper and in concomitant 
incorrect applications. These will be discussed in a sepa- 
rate subsection below. 



General comments on the relevance of AEs 

Before dwelling upon the more detailed discussion of 
the various models of modulated pulsation, we recall 
some of the relevant properties of the AEs. Most part of 
the subject discussed below can be summarized simply 



3 Stationary (constant amplitude) states correspond to steady multi- 
mode pulsations. These states are not interesting if only non-resonant 
radial modes are considered, because they are sparsely spaced and do 
not give rise to long-term modulation. 



as a warning against the straightforward application of 
the linear and nonlinear asymptotic pulsation results to 
the various pulsation states of the Blazhko modulation. 
Furthermore, the general framework of AEs enables us 
to test the ability of any new effect to cause modulated 
states. 

First of all, the concept of AEs is very general and 
can be applied whenever we would like to study the 
basic nonlinear properties of a dynamical system where 
the interaction goes on largely different time scales. For 
example, we can mention the long-term perturbations 
in celestial mechanics 13611 or processes on time scales 
considerably longer than the dynamical time scales in 
pulsating stars (such as mode switching or amplitude 
modulation, see JH, JH, |0]). Secondly, in many cases 
AEs are the only means to study complicated systems, 
where numerical simulations are either too expensive or 
impossible (e.g., nonlinear nonradial stellar pulsations). 

An important underlying assumption in deriving AEs 
for a given dynamical system, is that the nonlinearities 
are mild, meaning that they can be represented as a Tay- 
lor expansion in the neighborhood of the linear regime 
(that is usually rather easy to describe fairly accurately 
due to the separability of the time and spatial variations, 
and because of the linear superposition of eigenmodes 
in the general solution). If the above conditions are sat- 
isfied, then we start off by the following form of some 
physical quantity, e.g., of the displacement vector: 

8r(s,t) = fl (f)| (s)+fli(f)|'i(s) + ..., , (1) 

where s is the spatial coordinate, £, k is the linear non- 
adiabatic eigenvector of the normal-mode 'k' and a k (t) 
is the corresponding time-varying amplitude (complex, 
as all quantities in Eq. 1). The time-dependent functions 
are factorized to slow- and and short-varying parts and 
the AEs are derived for the slow-varying part. When 
Eq. (1) is substituted in the original partial differential 
equations that describe the fluid dynamics, we omit the 
short-oscillating terms and derive the AEs by demanding 
boundedness for the solution. 

The most important part of this process from the 
present context is that we start from the normal mode 
spectrum. This and the resonances will determine the ba- 
sic structure of the AEs representative for a given sys- 
tem. For example, if we consider the radial pulsation of a 
fully radiative stellar envelope in the gray atmosphere ap- 
proximation, we need to solve a system that is 3rd-order 
in time (two orders come from the time-dependence of 
the shell motion and another one comes from the en- 
tropy variation [see, e.g., fnj]). Therefore, this system 
leads to a 3 x iV-order eigenvalue problem for the nor- 
mal mode frequencies (N stands for the number of mass 
shells). Since the matrix of this eigenvalue problem con- 



stitutes of real elements, the eigenvalues are either pure 
real or come in complex conjugate pairs and a set of real 
ones. Due to the underlying physical setting, all radia- 
tive stellar models for classical pulsators display this lat- 
ter type of spectrum. The non-oscillating modes (due to 
their strong damping and different spatial variation) are 
usually considered to play no role in the behavior of the 
oscillating modes. The oscillating non-resonant normal 
mode spectrum leads to the following set of AEs in the 
lowest order approximation: 

dAi 

-± = A k (K k + Q k0 A 2 + Q kl Aj + ...) (2) 

= R k0 Al + R kl Aj + ... , (3) 

where, A k (t) and (p k (t ) stand for the slow-varying part of 
a k (t) in Eq. (1), i.e., the time-varying part of the mode 
behavior is described by A k (t)sm(co k t + (p(t)), with co k 
corresponding to the linear non-adiabatic eigenfrequency 
and K k to the linear growth rate. The coefficients {Q k j} 
and \R k i} are the non-resonant coupling coefficients. 
These are complicated expressions of the eigenfunctions 
and stellar structure, and rarely attempted to compute 
ab initio (see, however, B83I1 ). On the other hand, sim- 
ple physical considerations and numerical models show 
that each member of the set {Q k t} is negative. Similar 
statement may hold also for {R k i\ - i.e., Q9|] . 

It follows that for single-mode non-resonant pulsation, 
corresponding to the asymptotic state, after dying out 
the transients due to, e.g., the interaction with the other 
modes, the amplitude and the nonlinear period can be 
written in the following form: 

A\ = --^-+H.O.T. ( 4 ) 

Qkk 

2% 

P k (NL) = ^_^+H.O.T. , (5) 

a> k +R kk A 2 k 

where H.O.T. stands for 'higher order terms' (in ampli- 
tude). Since the coupling coefficients are slowly varying 
functions of the stellar parameters, the global properties 
of the single-mode non-resonant pulsation can be ap- 
proximated from the behavior of the linear growth rates. 

With these preliminaries, in the context of this review, 
we think that attempts trying to relate the observed tem- 
poral properties of Blazhko stars in terms of asymptotic 
and linear pulsation results are improper applications, 
because of the following two important reasons. 

■ Variations that take place in the large-amplitude 
asymptotic regime of stellar pulsation (i.e., in the 
observed variables) are related to the above (but not 
the same) types of equations and not to the ones 
described by the infinitesimal oscillations around 



the static state (i.e., the linear pulsational relations, 
e.g., van Albada and Baker jo t. 
• Equations (4) and (5) are irrelevant from the point 
of view of the Blazhko effect, since they are asymp- 
totic results whereas the time scale of the Blazhko 
effect is, in general, much shorter than the relax- 
ation time of the perturbations near the limit cy- 
cle s@ Furthermore, even if the time scales are com- 
parable, application of the asymptotic results (e.g., 
the nonlinear period change, see lf75ll . l34Tl ) at var- 
ious phases of the modulation, is in obvious con- 
tradiction with the original meaning of the above 
asymptotic values, since these kinds of applications 
assume instantaneous adjustment to the asymptotic 
state described by the non-modulated model. 

Based on the tight relation of the eigenspectrum of the 
linear problem and the type of AEs we can deal with 
in searching for systems exhibiting amplitude modula- 
tion, new/better physics can yield such a system, if they 
generate new modes that can interact with the radial pul- 
sation modes. In principle, convection might be a good 
candidate for generating such "new modes", but, so far, 
there are no reports on the instrumental nature of such 
modes, although current ID convection/pulsation mod- 
els do include several important properties of convection 
(H, El, JH, see also Buchler, these Proceedings). 
We do not see what kind of new treatment of convection 
can lead to a convective eigenmode spectrum that could 
give rise to AEs capable of leading to amplitude modula- 
tion (see some additional notes on the hypothetical role 
of resonant convective modes in j42tl ). 



The detailed modeling of the oblique rotator idea was 
largely motivated by the success of the same type of 
model by Kurtz [45] explaining the amplitude variation 
of the nonradial modes of the rapidly oscillating Ap stars. 
The first MORP model for the Blazhko effect by JMl 
was put in a more strict theoretical framework by 1 80] 
and 118111 as summarized in II64I1 . The essential result in 
these purely linear works is that under large-scale mag- 
netic fields all modes become nonradial. More specifi- 
cally, it was shown that in the case of a pure dipole mag- 
netic field the eigenmode corresponding to a radial mode 
in the non-magnetic case will be always deformed spa- 
tially by an / = 2 spherical harmonic component. It fol- 
lows that in the simplest setting we expect a quintuplet 
frequency structure around the radial pulsation compo- 
nent. The nice feature of this model is that the modula- 
tion components are scaled by the magnetic field strength 
and geometric factors (i.e., inclination angles) and these 
can be determined from the observed amplitude spec- 
trum. The two important obstacles that jeopardize this 
model are the following. 

• It follows from the above that the MORP model is 
purely geometric and the modulation is caused by 
rotation. Therefore we have strictly symmetric quin- 
tuplets both in amplitudes and in frequency separa- 
tions. There is no easy way to produce asymmetric 
modulation components in the frequency spectrajf] 

• The required strength of the dipole field is of the 
order of 1 kG. This is an order of magnitude larger 
that the one suggested by current measurements 
(see above). 



Magnetic oblique rotator/pulsator 

The first mentioning of an oblique-rotation mag- 
netic/pulsating configuration as a possible explanation of 
the the Blazhko phenomenon comes from yj] and j22tl . 
Although subsequent magnetic measurements by 10] and 
l59ll were controversial, the idea was further kept as a 
possible explanation due to some long-term modulation 
superposed on the 41 -day Blazhko cycle of RR Lyrae 
as being the manifestation of some magnetic star cycle, 
analogous to the sol ar cy cle |p78tl . Recent measurements 
of magnetic field by [ 17] strongly suggest that RR Lyrae 
itself should have a field lower than ~ 100 G. A simi- 
lar upper limit was obtained by 13911 by the spectropo- 
larimetric survey of 17 galactic field RR Lyrae stars that 
included at least 10 Blazhko stars. 



Nonradial resonant rotator/pulsator 

For an early, brief discussion and references on the 
possibility that nonradial modes might be instrumental 
in the Blazhko effect we refer to [73]. Starting in the 
early nineties, the idea has been revisited on a more 
fundamental basis by using updated linear numerical 
models and AEs (van Hoolst and coworkers in II 8411 and 
JUt], see also j42t] for a brief analysis of the relevant EAs 
and l28ll and ill 2ll for the discussion of AEs involving 
various nonradial modes). Further thorough elaboration 
of these ideas and specific extension to the Blazhko stars 
has been made by Nowakowski and Dziembowski in l54Tl 
(see also 12611 ). 

In the NRRP model we assume that the radial mode of 
interest (the FU or FO modes) are in 1 : 1 resonance with 
one of the not too strongly damped nonradial modes. 



4 The actual time scale depends on the type of perturbation and on the 
model (see, e.g., [ 76]), but for stable, dominantly single-mode pulsation 
in overtone 'k' this time scale is of the order of KZ 1 . 



5 In 1 64] there is a note on a possible break of this symmetry by an off 
axis dipole magnetic field, however, no detailed elaboration is given. 



Based on the models of l85ll and lf54Tl . we can strongly 
suspect that there are many nonradial modes of I = 1 that 
are in close resonance with the mode of interest and they 
are among the ones that are most viable for excitation 
(even though the code, on which this conjecture is based, 
does not include the effect of the interaction between 
pulsation and convection). The large number of available 
nonradial modes comes from the fact that in evolved stars 
the nonradial modes have high-order g-mode character in 
the deep interior. On the other hand, in the outer part of 
the stellar envelope they are similar to the p-modes and 
therefore, they have some chance to be spatially coupled 
with the radial modes (an important property when mode 
interaction is considered). 

From these encouraging linear properties it is sus- 
pected that the relevant AEs describing the resonant in- 
teraction between the radial and the nonradial mode(s) 
may lead to some stable state and thereby giving a chance 
for amplitude modulation either through some geometric 
or true dynamical interaction between the two classes of 
modes. In concentrating on the most viable interaction 
between the radial and the (7 = 1 , m = 0) or (I = 1 , m = 
±1) nonradial modes, it was shown by [54] that, in the 
asymptotic regime (i.e., for large time scales), the two 
types of interaction can be analyzed with the same set of 
equations valid for the 1 : 1 resonance with (I = 1 , m = 0) . 
Furthermore and very importantly, the most probable so- 
lution of the relevant AEs is the stationary one, with con- 
stant amplitudes and relative phase. This means that the 
interaction with the (7 = 1 ,m = 0) mode does not lead to 
amplitude modulation, because of the constant relative 
pha se. This results in the phase lock phenomenon (see 
112511 ). implying frequency synchronization, which, in the 
present case, means monoperiodic pulsation. Therefore, 
the NRRP model relevant for the Blazhko effect should 
consist of a 1:1 resonant stationary interaction with the 
(/ = l,m = ±1) modes and due to phase lock and stellar 
rotation, we observe a modulated pulsation. Although it 
is described briefly in ll54ll . and it is a simple derivation, 
we think that it is useful to see how the final frequency 
separation is related to the stellar rotation. 

Let the stationary nonlinear frequencies of the (/ = 
1 ,m = —1), (/ = 0) and (I = 1 ,m = +1) modes, re- 
spectively, in the co-rotating frame, be a)_, 0) and (0+. 
We adorn the corresponding linear values by zero su- 
perscript. Then, denoting the corresponding nonlinear 
changes by 8(0 , 8(0 and 8(0+, we have: 

(o = (o° + 8(0, (0+ = (0+ + 8(0± . (6) 

If the frequency distances between the central peak (cor- 
responding to the radial mode) and the side peaks (corre- 
sponding to the nonradial modes) are denoted by A±, we 
get: 

A_ = (o - (0- = (0° + 8(0 - (0° - 8(0 (7) 



A + = (o + - (o = a>l + 8(0+ -(0 -8(0 . (8) 

With the linear frequency mismatch of A© = 2(Oq — 
(0°_ — (0+, it follows that the difference between the two 
spacings can be expressed as: 

A+-A_ = -2A(0 + 28(0-8(0--8(0+ . (9) 

Then, using the stationarity condition and comparing 
the right-hand-side of the above expression with Eqs. 
(23) and (31) of fsill . we see that it is equal to zero, 
i.e., we have an equidistant triplet in the final state. 
Since the m = ±1 modes have the same radial structure, 
the corresponding nonlinear frequency changes are also 
the same, i.e., 8(0 = 8(0+. Therefore, the frequency 
separations can be expressed as follows: 

A_ = (o° + A(o) - (0 Q (10) 
A+ = (o+-(a° + A(o) . (11) 

Thus, the frequency separation A can be computed from 
the classical Ledoux formula, because: 

A=^± = ^^. (12) 
2 2 ' 

With Q. rotational angular velocity we have: (0+ — (0® ± 
CQ. + DD? + .... Since the nonradial modes are high- 
order g-modes, the coefficient C is equal to 0.5 in the 
asymptotic limit for / = 1 modes. The exciting outcome 
of all these is that if this model were viable, then the 
Blazhko period would be directly related to the rotational 
period of the star, so nonlinear period change would not 
influence the simple frequency splitting that we naively 
expect. In addition, the corresponding rotational rates are 
not in contradiction with the spectroscopically derived 
values [58], although for very short Blazhko periods, we 
might have problems [32]. There are also variables with 
multiple modulation periods 18711 . Obviously, the inter- 
pretation of these objects in terms of rotation modulation 
is even more difficult. 

As far as the amplitudes of the modulation compo- 
nents are concerned, they will be symmetric both for 
physical and for geometrical reasons (same type of spa- 
tial dependence and rotation). As a result, the main rea- 
son of the fall of the NRRP model can be summarized as 
follows. 

■ In the basic model described above the modulation 
components have strictly equal amplitudes. This 
property is in clear contradiction with the observa- 
tions that show the overwhelmingly opposite behav- 
ior of Blazhko stars. Perhaps the consideration of 
additional nonradial modes could lead to asymmet- 
ric amplitudes, but this would require the introduc- 
tion of additional unknown parameters and the re- 
sult would not be verifiable. 



However, there seems to be also a deeper problem with 
our current understanding of nonradial mode coupling 
in evolved stars. In ifsll the question has been posed if 
the condition of mild nonlinearity (one of the basic as- 
sumptions behind the applicability of the AE formalism) 
is satisfied for both the radial and the horizontal com- 
ponents of the nonradial modes. The most essential part 
of the reasoning is that the horizontal component of the 
displacement vector is proportional to the Brunt-Vaisala 
frequency (see their Eqs. (1), (2)). This latter quantity is 
a very sensitive function of the internal structure and can 
reach very high values near the core of evolved objects. 
Figure 1 shows the variation of the relative Brunt-Vaisala 
frequency as a function of the distance from the star's 
center. We see the large difference between (near) Main 
Sequence and Horizontal Branch stars. In the latter the 
Brunt-Vaisala frequency can be higher by nearly two or- 
ders of magnitude in the deep interior than in MS stars. 
The large value implies also large horizontal displace- 
ment (several hundred times larger than that of the radial 
component). Based on this result, the authors also make 
an estimate on the expected size of the critical surface 
amplitudes of the nonradial modes, where nonlinearity 
becomes strong (i.e., advective terms start to dominate 
the fluid motion). For / = 1 modes in standard RR Lyrae 
models they obtained an amplitude of 0.02 mag. This is 
considerably smaller than the generally observed modu- 
lation side lobes in the frequency spectra of the Blazhko 
stars. These results suggest that the standard AE formal- 
ism that assumes mild nonlinearity might not be appli- 
cable for studying nonradial mode interaction in evolved 
stars. 
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FIGURE 1. Variation of the Brunt-Vaisala frequency for dif- 
ferent stellar models as a function of the stellar radius (adapted 
from arXiv: astro-ph/030 1 504V 1 , see also (H]). 



Dynamo+convection+pulsation + criticism 

The purpose of this subsection is to draw attention to 
the inconsistencies occurring in a recent idea of Blazhko 
effect presented in I75I1 . Some of the problems concern- 
ing this idea have already been mentioned in earlier sec- 
tions. Here we give a brief, but more complete review. 

• First of all, the idea lacks the level of elaboration 
observed in all other often cited models (MORR 
NRRP and 2:1 resonance models - see next sec- 
tion). There is no prediction on the expected mod- 
ulation period or the level of modulation. No men- 
tioning is made on the generic case of unequal mod- 
ulation side lobes (one serious problem that makes - 
among others - the MORP and NRRP models fail). 

• Secondly, there is a misinterpretation of the non- 
linear hydrodynamical results concerning the rele- 
vance of the asymptotic periods in explaining the 
temporal periods appearing in one type of mathe- 
matical forms of the Blazhko modulation. This topic 
is discussed already earlier, so we just mention that: 
(a) All nonlinear radially pulsating stellar models 
show small period increase in their final single- 
mode states (except for @], but see JiH], HU and 
also Smolec, private communication); (b) These pe- 
riod changes have nothing to do with the temporal 
period (as defined by the time derivative of the total 
phase) of the Blazhko modulation. Nonlinear period 
changes are asymptotic and reaching these states 
takes a long time; (c) The amplitude-dependence of 
the nonlinear periods as predicted by the AEs has 
also an asymptotic meaning, and cannot be applied 
to explain the possible correlation between the tem- 
poral phase (i.e., period) and amplitude values of 
the Blazhko modulation. 

• Thirdly, the heart of the idea is an assumed peri- 
odic build up and destruction of a magnetic field 
via the interaction with the turbulent convection. 
No mentioning is made on the expected strength 
of this field and why this would lead to an ampli- 
tude modulation. The way it is presented, sounds 
like if the assumption was that pulsation is work- 
ing in an environment determined by the interaction 
between convection and turbulence and there is no 
feedback between them. If this is so, then the idea 
lacks the important effect of convection-pulsation 
interaction, and, due to the assumed magnetic field- 
turbulence interaction, the pulsation-magnetic field 
feedback is also missing. In very general terms, it is 
unclear why this highly complicated system would 
lead to AEs that are capable of giving rise to modu- 
lated pulsation. (We recall that the eigenmode spec- 
trum of a dynamical system will determine the pos- 
sible forms of AEs, so not knowing the linear prop- 



erties of the full linear system, we cannot predict 
what happens in the nonlinear regime). 

In conclusion we think that the complex idea depicted 
in [ 75] could be verified only if some sort of deeper phys- 
ical modeling were possible, first in the regime of lin- 
ear magnetohydrodynamical convective/pulsation mod- 
els. This seems to be quite a complicated but perhaps not 
an impossible task. 

A RESONANT TOY MODEL - FOR 
ENTERTAINMENT 

This section is by no means to suggest that the model 
presented could be seriously considered as one of the 
possible explanations of the Blazhko effect. Our goal 
is just to introduce a non-conventional hydrodynamical 
model that shows an amplitude- and phase-modulation, 
resembling to what we see in the observations. Since this 
is the first time such a toy model is presented, we think it 
may be interesting and will perhaps encourage others to 
pursue some other types of non-standard models. 

In 1986 Moskalik [48] suggested that a 2:1 resonance 
between the fundamental and third overtone modes 
in RR Lyrae stars might cause real dynamical ampli- 
tude/phase modulation. The idea follows a model of 
double-mode pulsation 1I25I1 which was thought to be en- 
couraging at that time, and which utilizes the same type 
of resonance. The equations governing the mode interac- 
tion under this resonance are as follows: 

— - = Ai(Ki+qnA\ + q n Al + quAl) 
at 

+ C1A1A3 cos4> (13) 

= A 2 (K 2 + q 2l Aj + q 2 2Al + q 2i Al) (14) 

at 

— = A 3 (K3 + q 3l Ai + q 32 A^ + q 33 Ai) 

- C 3 A\ cos4> (15) 

- = ri A\ + r3 Al 

+ Aa)+r^--2CiA 3 Jsin<5 . (16) 

Here we omitted non-adiabaticity in the resonant cou- 
pling Cfc. In the above notation < and Q > for all 
(i,j,k). From the increase of the nonlinear periods in the 
case of non-resonant hydrodynamical models it follows 
that the parameters are also negative. The frequency 
mismatch is computed from the linear eigenfrequencies 
as Ao) = C0 3 -2<0\. 

The most general stationary solution of these equa- 
tions is the so-called "bump" solution with A 1 ^ 0, A 2 = 



and A3 7^ 0. In this case we have a single-period pul- 
sation due to phase lock and the higher frequency over- 
tone causes a bump appearing, e.g., on the radial velocity 
curve. This solution is thought to be the underlying phys- 
ical model of the bump Cepheids (Simon and Schmidt, 
Buchler and coworkers, see JH, M, iflHl ). 

Yet another stationary solution of the above equa- 
tions was proposed by Dziembowski and Kovacs 12511 
to model steady double-mode pulsation. In this model 
the resonance between mode '1' and '3' leads to the de- 
crease of the amplitude of mode ' 1 ' relative to its single- 
mode non-resonant value. This leads to a lower stability 
of the above 'bump' solution against perturbations to- 
ward mode '2'. If the resonance is strong enough (Acq 
is small) and some other, not too restrictive conditions 
are also satisfied, mode '2' will grow and this will even- 
tually lead to a steady triple-mode pulsation (with two 
'essential' periodicity in the system, since mode ' 1 ' and 
'3' remain in resonant frequency synchronization). Al- 
though the mechanism was shown working also in direct 
hydrodynamical simulations iPTill . the required resonance 
condition is not satisfied in standard RR Lyrae modelsQ 

In the work of Moskalik l48tl two additional impor- 
tant simplifications were made in the 2:1 resonance sys- 
tem given by Eqs. (13)— (16): (a) the interaction with the 
non-resonant FO mode (Eq. 14) was omitted; (b) sat- 
uration in the third mode was discarded (i.e., all {^3,} 
were set equal to zero), but it was assumed that the mode 
was linearly damped (k 3 < 0). With the simplified sys- 
tem the condition of stability of the resonant solution 
was addressed by allowing wide ranges of parameters, 
more importantly those of the frequency mismatch Ao 
and the damping rate of the resonant overtone K 3 . The 
parameter survey disclosed a wide range of possible lim- 
iting states. For high damping rates of mode '3' the sta- 
ble two-mode fixed point solution was the only solution. 
The situation became rather complex when the damping 
rate was lowered (see his Fig. 1). In general, for a very 
a mild damping of mode '3', no stable fixed point was 
found for very wide a range of Acq. Since, for standard 
RR Lyrae models, the resonant mode is the 3rd over- 
tone and it is not too mildly damped, Moskalik investi- 
gated this regime in more detail. Here, the modulated so- 
lution was achieved for non-exact resonances (Ao 7^ 0) 
with some additional conditions put on the non-resonant 
period change terms (parameters {r^} in Eq. 16). Un- 
fortunately, all these seem to be rather restrictive, since 



6 In the case of Cepheids we do have such a resonance but it occurres 
in the bump regime - near 10 days period - that is not where most 
of the beat (double-mode) Cepheids can be found. Nevertheless, the 
destabilizing effect of the 2:1 resonance in respect of the single-period 
bump solution was also demonstrated for this class of models by 
hydrodinamical simulations by fTTIl . l49ll . 



TABLE 3. Stellar and model construction parameters of the Blazhko toy model. 



Code: 
Zoning: 

Stellar parameters: 
Composition: 



Castor/Stellingwerf radiative, LANL opacities 

^sheii = 100, Mnz = 30* 

M/Mq = 0.90, L/Lq = 40.0, r eff = 6400 K 

X = 0.760, Z = 0.001 throughout, except in 1.0£5 < T < 1.0£6, where 
X = 0.661, Y = 0.239, Z = 0.100, with ramped Z/opacity-enhancement of 
100 and 20, respectively. 



Numerical viscosity: CQ = 4.0, a = 0.005 



Number of mass shells from the surface down to the Hydrogen ionization zone at T = 1 1000 K. 



TABLE 4. LNA properties of the Blazhko toy model. 



Type 


Po 


Pi 


Pi 


P3 


P2/P0 




no 


m 


m 


H3 


Pi /Pi) 


NOT ENH*_ 


0.42 


0.32 


0.26 


0.21 


0.62 




0.01 


0.03 


0.03 


0.00 


0.76 


ENHj. 


0.53 


0.33 


0.26 


0.21 


0.49 


0.02 


0.03 


0.04 


-0.01 


0.62 



* Standard linear models (no opacity/metal enhancements). 
^ Opacity/metal enhanced models as described in the text. 



the large number of various hydrodynamical simulations 
performed during the past 23 years did not show signs of 
modulated variations. 

In contemplating the disappearance of stable fixed 
point solutions in the mildly damped case, we decided to 
test the case when this happens in real hydrodynamical 
models. We searched for a means to bring one of the 
overtones in 2:1 resonance with the fundamental mode 
and make it less damped (or if possible, even excited, 
since the more physical hydrodynamical models do have 
the saturation terms {qn} that have been neglected in 
l48tl and they will always increase the effective damping 
in mode '3'). 

To create RR Lyrae models with the above proper- 
ties, we tried the 'good old-fashioned' opacity/heavy el- 
ement enhancement l66ll . [88] [18]. The required size of 
enhancement is rather large, so there is no good physi- 
cal reason to introduce such change in the models [j Our 
rather vague excuse for doing this is the assumption that 
there was an element enhancement in the earlier phase of 
evolution (e.g., during the ascent to the first giant branch, 
due to the capture of a companion of high heavy ele- 
ment content - see JHH). Somehow mixing was avoided 
and most of the heavy elements settled in the deep inte- 
rior (below the outer convection zone). As a result, they 
do not affect the measured iron content for these typi- 



cally low-metal stars Admittedly, this whole scenario is 
highly unlikely, but the present toy model is not intended 
to be too physical. 

For computing the purely radiative static, linear non- 
adiabatic (LNA) and nonlinear models, we used the same 
code as in our earlier works [ 10]. For a comfortable scal- 
ing of the opacity, we used Stellingwerf's ll72ll analytical 
formula. For an easier overview, in Table 3 we summa- 
rized the stellar and model construction parameters. We 
note that the large opacity enhancement is relative to the 
Los Alamos opacities, prior to the OPAL/OP revision of 
these l6lll . 16211 . rt30ll . Therefore, in respect of the, e.g., 
OPAL opacities the increase is 'only' a factor of five or 
so. In addition to the model presented we computed sev- 
eral/many others and the properties to be described be- 
low can be regarded as generic. 

The LNA properties of the model are shown in Ta- 
ble 4. For comparison, we also give the corresponding 
non-enhanced properties. We see that the largest changes 
due to Z/OP enhancement occurred in the FU mode. The 
drastic change in the period has resulted the reach of 
near 2:1 resonance condition between Po and P2 with a 
moderate increase in the growth rate of the second over- 
tone. These are very favorable changes for the antici- 
pated condition of amplitude modulation within the res- 
onant model of [48] (not so for the standard horizontal 
branch models, since the values of P\/Pq and Pq are out 
of range for the double-mode stars and for the PLC rela- 
tions ll43Tl ). 

We initialized the static model with a velocity pertur- 
bation corresponding to the fundamental mode and in- 
tegrated in time until the asymptotic state was reached. 
Figure 2 shows the frequency spectra of the bolomet- 
ric light variation of the last 500 Po cycles. Succes- 
sive prewhitening of the time series revealed the vari- 
ous combination frequencies of Vo = 1.88939 d _1 and 
v,„ = 0.04684 d , corresponding to the fundamental 
mode and to the modulation frequencies, respectively. It 



7 See however the 'Iron Project' [ 56] and l47h for some renewed efforts 
on the revisitation of the opacities of complex heavy elements. 



8 It is interesting to note that in the Blazhko star TY Gruis there are 
signatures of very large enhancements in some heavy elements (see 
1 60 ]), but this is most probably some peculiarity of this object rather 
than a general feature of the Blazhko stars. 
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FIGURE 3. Bolometric light curve of the resonant toy 
model. Panels from top to bottom show the light curve folded 
by the pulsation period, close-up of the time series and longer 
stretch of it, covering three modulation periods. 



FIGURE 2. The frequency spectra of the bolometric light 
curve of the resonant toy model (see Tables 3 and 4). Each 
spectrum is normalized to the highest peak. Numbers in the 
insets of the panels show: the prewhitening order, the peak 
frequency (in [d ]), the corresponding amplitude (in [mag]) 
and the signal-to-noise ratio (which is given as 99.9, if it is 
greater than 100). 



is important to note that the components are in very tight 
numerical relations. Alternatively, and perhaps more nat- 
urally, we can represent all frequencies as linear combi- 
nations of Vo above and V2 = 3.82559 d _1 , correspond- 
ing to the second overtone. This representation is very 
similar to the one advocated by [|5|] and ((J. Indeed, if 
the asymptotic state is an FU/SO double-mode pulsa- 
tion rather than a genuine amplitude modulation, then the 
higher frequency modulation component at Vo is V2 — Vo, 
whereas the much lower amplitude one at the lower fre- 
quency side corresponds to 3Vo — V2, in line with the ex- 
pected lower amplitude for such a high order combina- 
tion. 

The bolometric light variation is shown in Fig. 3. In 
the middle panel we see the movement of the bump due 
to the presence of the second overtone. The folded light 
curve exhibits the combination of phase and amplitude 
variations. All these are also observable in real Blazhko 
stars. 

We think that the the above model works the way a 
toy model should; it is certainly entertaining, but requires 
most probably unacceptable torture of the stellar struc- 



ture. Further temptations for additional contemplation of 
a possible viability of the idea is hampered also by the 
lack of models with nearly symmetric modulation com- 
ponents. 

CONCLUSION AND PROSPECTS 

In spite of the renewed effort during the past 15 years, 
so far none of the observational and theoretical studies 
led to a breakthrough in our lack of understanding the 
Blazhko phenomenon. More than a thousand RR Lyrae 
stars are known in various stellar systems that show peri- 
odic amplitude- and phase-modulations. Importantly, the 
phenomenon has been detected recently also in Cepheids 
i52ll with a seemingly quite high incidence rate (at least 
for some classes of Cepheids). There are observational 
projects (Jurcsik and coworkers, Kolenberg and cowork- 
ers) that study more individual cases and see if there is 
anything in the light and color variations that might give 
us a clue on the underlying mechanism. Sparsely though, 
but very importantly, there are also spectroscopic works 
aimed at the detection of magnetic field in RR Lyrae stars 
01, lUl. Theoretical works JH, JH that targeted the 
two main (magnetic oblique rotator/pulsator and nonra- 
dial resonant rotator/pulsator) models have failed mostly 
for being unable to interpret the unequal modulation 
components in the frequency spectra and also for other 
observational and theoretical reasons (strong nonlinear- 
ity of the nonradial modes and apparent absence of 



strong magnetic field). 

It seems to be that the inspiration for new types of 
models should come from the observational side. Al- 
though current ground-based observations can also tar- 
get searches for additional frequency components in 
the light variation in the millimagnitude level, the real 
breakthrough will come from the current space projects 
CoRoT and Kepler. First signs of this event have already 
been witnessed from the preliminary results coming from 
CoRoT (Chadid et al., these Proceedings). Whether the 
appearance of additional small amplitude components 
will be a common feature of all Blazhko stars or this 
feature remains to be a rare occurrence, is still to be 
seen. Obviously, the presence of low-amplitude modu- 
lated group of frequencies would be suggestive of a so 
far not investigated effect of the dynamical interaction 
of several/many nonradial and radial modes 112711 . It is 
also interesting how the incidence rate of the Blazhko 
phenomenon will change as the signal-to-noise ratio of 
the measurements increases. There are already signs that 
the phenomenon is much more common than we thought 
(Jurcsik, these Proceedings). Unfortunately, the number 
of the spectroscopic works is far below that of the pho- 
tometric ones. We think that line profile analyses (even 
if they were based on snapshot observations) would be 
very important for future modeling. An upper limit on 
the size of the nonradial component of the velocity field 
would be instrumental in considering nonradial modes 
either as members of a periodic energy sink, or retain 
their leading role in establishing the amplitude modula- 
tion, or reject the nonradial mode interaction altogether 
and search for other mechanisms. 
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